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Abstract. We continue to study a model of disordered interface growth in two dimensions. The 
interface is given by a height function on the sites of the one-dimensional integer lattice and 
grows in discrete time: (1) the height above the site x adopts the height above the site to its left 
if the latter height is larger, (2) otherwise, the height above x increases by 1 with probability 
Px- We assume that are chosen independently at random with a common distribution F, and 
that the initial state is such that the origin is far above the other sites. Provided that the tails 
of the distribution F at its right edge are sufficiently thin, there exists a nontrivial composite 
regime in which the fluctuations of this interface are governed by extremal statistics of p^- In 
the quenched case, the said fluctuations are asymptotically normal, while in the annealed case 
they satisfy the appropriate extremal limit law. 
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FLUCTUATIONS IN THE COMPOSITE REGIME 
OF A DISORDERED GROWTH MODEL 

Janko Gravner, Craig A. Tracy, Harold Widom 

1. Introduction. 

Disordered systems, which are, especially in the context of magnetic materials, often referred 
to as spin glasses, have been the subject of much research since the pioneering work in the 1970s. 
The vast majority of this work is nonrigorous, based on simulations and techniques for which a 
proper mathematical foundation is yet to be developed. (See [MPV] for early developments and 
[Tal] for a nice overview of the mean field approach.) As a result, there is a large number of new 
and intriguing phenomena observed in these models which await rigorous treatment. Among the 
most fundamental of issues are the existence and the nature of a phase transition into a glassy 
or composite phase: below a critical temperature, the dynamics of a strongly disordered system 
becomes extremely slow with strong correlations, aging and localization effects and possibly 
many local equilibria. We refer the reader to [NSv] and [BCKM] and other papers in the same 
volume for reviews and pointers to the voluminous literature and to [NStl] and [NSt2] for some 
recent rigorous results. In view of the difficulties associated with a detailed understanding 
of realistic spinglass systems, other disordered models have been introduced, which are more 
amenable to existing probabilistic methods. 

One of the most successful of such (deceptively) simple models is the one-dimensional random 
walk with random rates [FINl]. In this model, the walker waits at a site x G Z for an exponential 
time with mean Tx before jumping to either of its two neighbors with equal probability. The 
disorder variables are i.i.d. and quenched, that is, chosen at the beginning. Provided that the 
distribution of has sufficiently fat tails, namely, if P{tx > t) decays for large t as as with 
a < 1, the walk exhibits aging and localization effects ([FINl], [FIN2]). Various one-dimensional 
voter models and stochastic Ising models at zero temperature can be explicitly represented with 
random walks. This connection has been explored to demonstrate glassy phenomena such as 
aging and chaotic time dependence ([FINl], [FINS]). The positive temperature versions of such 
results remain open problems, even in one dimension. 

In contrast with models which are exactly solvable in terms of random walks and are by now 
a classical subject in spatial processes ([Gril], [Lig]), techniques based on the RSK algorithm 
and random matrix theory have entered into the study of growth processes only recently ([BDJ], 
[Johl], [Joh2], [BR], [PS], [GTWl]). The purpose of this paper is to employ these new methods 
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to prove the existence of a pure phase and a composite phase in a disordered growth model. 
It has been observed before in similar models [SK] that the role of temperature is for fiat 
interfaces apparently played by their slope. In our case, the initial set is very far from fiat and 

"temperature" is measured instead by the macroscopic direction (from the origin) of points on 
the boundary. We identify precisely the critical direction and demonstrate that the fluctuations 
asymptotics provide an order parameter that distinguishes the two phases. We emphasize that 
a hydrodynamic quantity, the asymptotic shape, has a discontinuity of the first derivative at the 
transition point, at which the shape changes from curved to flat. However, this docs not signify 
the existence of a new phase as kinks are common in many random growth models [GG], thus 
a finer resolution is necessary. 

The particular model we investigate is Oriented Digital Boiling (ODB) (Feb. 12, 1996, Recipe 
at [Gri2], [Gra], [GTWl], [GTW2]), arguably the simplest interacting model for a growing 

interface in the two-dimensional lattice Z^. The occupied set, which changes in discrete time 
i = 0, 1, 2, . . . , is given by At = {{x, y) : x ^ 7i,y < ht{x)}. The initial state is a long stalk at 
the origin: 



while the time evolution of the height function ht is determined thus: 

ht+i{x) = max{ht{x - l),ht{x) + £x,t}- 

Here e^^t s^re independent Bernoulli random variables, with P{£x,t = 1) = Px- Although this 
model is simplistic, note that it does involve the roughening noise (random increases) as well 
as the smoothing surface tension effect (neighbor interaction), the basic characteristics of many 
growth and deposition processes. (See Sections 5.1, 5.2 and 5.4 of [Mea] for an overview of 
simple models of ODB type as well as some other disordered growth processes.) 

We will assume, throughout this paper, that the disorder variables Px are initially chosen at 
random, independently with a common distribution F{s) = P{px < s). We use ( ■ ) to denote 
integration with respect to dF and label by p a generic random variable with distribution F. 

It quickly turns out ([GTWl]), that fluctuation in ODB can be studied via equivalent increas- 
ing path problems. Start by constructing a random m x n matrix A = A(F), with independent 
Bernoulli entries Si^j and such that P{£ij = 1) = Pj, where, again, pj = p are i.i.d. Label 
columns as usual, but rows started at the bottom. We call a sequence of I's in A whose posi- 
tions have column index nondecreasing and row index strictly increasing an increasing path in 
A, and denote hy H = H{m, n) the length of the longest increasing path. Then, under a simple 
coupling, ht{x) = H(t — x,x + 1) ([GTWl]). Thus we will concentrate our attention on the 
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random matrix A rather than the associated growth model. Prom now on we will also replace 
Pi with its ordered sample, so that Pi > P2 > ■ ■ ■ > Pn (see section 2.2 of [GTWl]). 

We initiated the study of ODB in a random environment in an earlier paper ([GTW2]), from 
which we now summarize the notation and the main results. Throughout, we denote by b the 
right edge of the support of dF and assume it is below 1, i.e., 

b = mm{s : F{s) = !}<!. 

Moreover, we fix an a > and assume that n = am. (Actually, n = [am\, but we omit 
the obvious integer parts.) As mentioned above, we can expect different behaviors for different 
slopes on the boundary of the asymptotic shape, which translates to different a's. To be more 
precise, we define the following critical values 

\ -1 

. P 

an 



a' 



\ ~ p 
p{l-p) 



{b-pf 

Note that the second critical value is nontrivial, i.e., a'^ > 0, iff — p)~^) < 00. Next, define 
c = c(a, F) to be the time constant, c = c{a, F) = lim H/m, which determines the limiting 

m — *oo 

shape of At, namely limAt/t, as t ^ 00. In Theorem 1 of [GTW2], it was found that c exists 
a.s. and is given by 

' b + a{l-b){p/{b-p)) , iia<a'^, 
c{a, F) = < a + a(l — a) {p/ (a — p)) , if aj, < a < Oc, 
1, if CKc < a. 

Here a = a{a, F) G [b, 1] is the unique solution to a (p(l —p)/{a — p)"^) = 1. 



In [GTW2], we also determined fluctuations in the pure regime < a < Oc- (The deter- 
ministic regime ac < a has no fluctuations.) The annealed fluctuations ([GTW2], Theorem 2) 
about the deterministic shape c grow as ^Jm and are asymptotically normal: 

H — cm d 



iV(0,l) 



as m — 00, where Tq = Var((l — a)p/ (o — p)). 



By contrast, quenched fluctuations conditioned on the state of the environment grow more 
slowly, as m^^^, and satisfy the F2-distribution known from random matrices ([TWl], [TW2]). 
To formulate this result, we let rj = pj/{l — pj), define u„ to be the solution of 



(1.1) ^f^ "1 = ^ 

^ ^ n ^ (1 + rjuy {u - 1)2 
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which Ues in in the interval {—r-^ "^,0). This solution exists provided that an ^J2^=i''~j < 1 
which holds a.s. for large n as soon as a < ac- Next, set c„ = c{un) where 



1 " 

(1.2) c{u) = rr— --Yt 

^ ' ^ ^ 1 - w n ^ 1 




Then ([GTW2], Theorem 3) there exists a constant 7^ so that 



P 



H — c„m 





as m — > oo, almost surely, for any fixed s. 

For fluctuation results in this paper we need to impose some additional assumption on F, 
which are best expressed in terms of G{x) = 1 — F({b — x)—), the distribution function for b — p. 
First we list our weaker conditions: 

(a) li x,y ^ and x ^ y, then G{x) ~ G{y). 

(b) If ^ and X = then G{x) = 0{G{y)). 

(c) As X ^ 0, G{x) = o{x^/\ogx~^). 

Our stronger assumptions on F require that there exists a 7 > so that: 
(a') The function G{x)/x^ is nonincreasing in a neighborhood of x = 0. 
(b') G{x) = 0(x^/log'' x~^) as X ^ for some > 27 + 4. 

If a'^ > 0, then automatically G{x) = o(x^) as x — 0. The stronger assumptions thus do 
not require much more: for nicely behaved G they amount to G{x) = 0(x^/ log"^ x~^) for some 
v > 8. The quenched and annealed fluctuations are now determined by the next two theorems. 

Theorem 1. Assume that < a < a'^, let 



and let $ be the standard normal distribution function. If (a)-(c) hold, then for any fixed s, as 
m — ^ 00, 



Here, the convergence is in probability if (a)-(c) hold, and almost sure if (d ) and (b'J hold. 





5 



Theorem 2. Assume that < a < a^, and that (a)-(c) hold. Then, for any fixed s 

P [H < cm — {I — a/a'^)mG~^{s/n) \ pi,...,Pn) e~* 
in probability. In particular, 

P{H <cm-{l-a/a'^)mG-^{s/n)) ^ e''. 

Throughout, we follow the usual convention in defining G~^{x) = sup{y : G{y) < x} to be 
the left continuous inverse of G, although any other inverse works as well. 

Assume, for simplicity, that, as x — 0, G{x) behaves as for some rj > 2. Then, in contrast 
with the pure regime, the annealed fluctuations in composite regime scale as m^~^^'^, while the 
quenched ones scale as m^/^. In fact, this can be guessed from [GTW2]. Namely, as explained 
in Section 2 of that paper, the maximal increasing path has a nearly vertical segment of length 
asymptotic to (1 — a/a'^)m in (or near) the column of A which uses the largest probability pi. 
Therefore, this vertical part of the path dominates the fluctuations, as the rest presumably has 
o{y/rn) fluctuations. (These are most likely not of the order exactly m^/'^ as they correspond 
to the critical case a = a'^. The precise nature of the critical fluctuations is an interesting open 
problem.) The variables in the pi-column are Bernoulli with variances about 6(1 — 6), thus the 
contribution of the vertical part to the variance is about (6(1 — 6)(1 — a/a^)m)^/^ = r^/n. The 
annealed case then simply picks up the variation in the extremal statistic pi . 

Simple as the above intuition may be. Theorems 1 and 2 are not so easy to prove and 
require considerable additional technical details. We also note the mysterious correction 2T^fn 
in Theorem 1 for which we have no intuitive explanation. 

The fluctuations results in [GTW2] and the present paper thus sharply distinguish between 
two different phases of one particular growth model. Nevertheless, it seems natural to speculate 
that this phenomenon is universal in the sense that it occurs in other one-dimensional finite 
range dynamics of ODB type, started from a variety of initial states. Indeed, such universality 
has been established in other random matrix contexts [Sos] . Fluctuations of higher-dimensional 
versions seem much more elusive; it appears that a glassy transition should take place, but the 
fluctuation scalings could be completely different. 

To elucidate, we present some simulation results. In all of them, we start from the flat 
substrate /iq = and use -F(s) = 1 — (1 — 2s)'', so that 6 = 1/2. It is expected that, as t] 
increases, the quenched fluctuation experience a sudden jump from 1/3 to 1/2. We simulate two 
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dynamics, the ODB and the two-sided digital boihng (abbreviated simply as DB), given by 

ht+i{x) = max{ht{x - l),ht{x + l),ht{x) + e^^J. 

The top of Figure 1 illustrates the ODB on 600 sites (with periodic boundary), run until time 
600. The occupied sites are periodically colored so that the sites which become occupied at the 
same time are given the same color. On the left, r] = 1 (i.e., p is uniform on [0, 1/2] and = 0), 
while 7] = 3 (and > 0) on the right. The darkly colored sites thus give the height of the 
surface at different times and provide a glimpse of its evolution. In the pure regime {rj = 1), the 
boundary of the growing set reaches a local equilibrium ([SK], [BFL]), while in the composite 
regime {rj = 3) the boundary apparently divides into domains, which are populated by different 
equilibria and grow sublinearly. This is the mechanism that causes increasing fluctuations. The 
bottom of Figure 1 confirms this observation; it features a log-log plot of quenched standard 
deviation (estimated over 1000 independent trials) of /it(0) vs. t up to t = 10 000. The rj = 1 
case is drawn with +'s and the ?? = 3 case with x's; the two least squares approximations lines 
(with slopes 0.339 and 0.517, respectively) are also drawn. We note that the asymptotic speed 
of this flat interface is known: limt^oo ht{0)/t = sup^^ol*^ + l)c(a). Here is the reason: if ODB 
dynamics h], ht start from initial states Hq, = supj /iq, respectively, and are coupled by using 
the same coin flips ex,t, then ht = supj hi for every t. 




Figure 1. Evolution and quenched deviation in the two phases of disordered ODB. 
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Perhaps surprisingly, it appears that the phase transition in the DB does not occur at 77 = 2, 
and in general the delineation is much murkier. At this point, we cannot even eliminate the 
possibility of continuous dependence of fluctuation exponent on rj. In Figure 2, we present the 
results of simulations for r] = 0.2 (left) and r] = 1 (right). The top figures only show evolution 
near time t = 5000, as no difference is readily apparent at earlier times. The plot of quenched 
deviations is analogous to the one in Figure 1, with the least squares slopes 0.395 {r] = 0.2) and 
0.49 (t/ = 1). 




0123456789 10 

Figure 2. Evolution and quenched deviation in disordered DB. 

The organization of the rest of the paper is as follows. Section 2 reviews the set-up from 
[GTWl, GTW2], in Section 3 we prove the relevant asymptotic properties of the order statistics 
and of the solutions of (1.1) and (1.2), and demonstrate how Theorem 2 follows from Theorem 
1. Section 4 is a detailed analysis of the asymptotic behavior of steepest descent curves. The 
proof of convergence in probability in Theorem 1 is then concluded in Section 4. Finally, Sec- 
tion 5 strengthens the results of Section 3 (under the stronger conditions) so that almost sure 
convergence is implied. 
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2. The basic set-up. 

We recall how we approached these problems in [GTW1,GTW2]. The starting point is the 
identity 

P{H<h) = det {I -Kh), 
where is the infinite matrix acting on ^^(Z"*") with (j, fe)-entry 

oo 

e=o 

The subscripts denote Fourier coefficients and the functions (p± are given by 

n 

= + rjz), ^.(z) = (1 - z-')-^. 
The matrix is the product of two matrices, with {j, A;)-entries given by 



1 r " 

= — / Hil + r.z) {z-irz-^+^+^+^dz, 

1 r 

{^./^+)n+j+k+i = ^ / Y[{l + r,z)-^ {z-iy-^z^-^-^-^-^dz. 

The contours for both integrals go around the origin once counterclockwise; in the second integral 
1 is on the inside and all the —tJ^ are on the outside. 

li h = Cnm + h' WQ have 

(2.1) {^+/v.).h-j-k-i = ^ / V'(^) z'^'^^^^ dz, 

(2.2) {<p./^+),+j+k+i = ^.J i^iz)-' z-^'-^-'^-^ dz, 
where 

n 

■il,{z) = J[{l + rjz) (z-l)"^z-(i-'=")'". 

.7 = 1 

The idea is to apply steepest descent to the above integrals. If a{z) = log ip{z) then 

(2.3) a'{z) = ^ V + -^ + ^ 

^ ' ^ ' n 1 + r^z z-l z 

j=i ■> 

and, with n„ and c„ as defined above, cj'{un) = cr"{un) = 0. The steepest descent curves both 
pass through n„. As n — > oo the zeros/poles —tJ^ accumulate on the half-line {-co, ^] where 
^ = 1 — . In the pure regime the points Un and the curves are bounded away from this 
half-line, behave regularly and have nice limits. However in the composite regime the points 
and curves come very close to 4, their behavior is not so simple, and we apply steepest descent 
not quite as described. 
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3. Preliminary lemmas I: properties of Pn, Un, and c„. 

Until Section 5, wc assume that all limits are in probability, unless otherwise indicated. To 
prove the first part of Theorem 1 and Theorem 2, we thus assume that (a)-(c) hold. 

We let Qj = b — pj , so that qi, ■ ■ ■ ,qn are chosen independently according to the distribution 
function G, then ordered so that qi < q2 < ■ ■ ■ < Qn- 

Let ti < t2 < ■ ■ ■ < tn he an ordered sample of i.i.d. uniform (0, 1) random variables. Then 
we may construct the G-sample by setting qj = G~^{tj). We will also use the well-known fact 
that, given tj, the conditional distribution of ti, . . - tj-i is that of an ordered sample of j — 1 
uniforms on [0,tj]. 

Lemma 3.1. There exist a positive constant Ci so that x < G{G~^{x)) < x/ci for x G (0, 1). 
Moreover, G{G~^ (x)) ~ x as x — ^ 0. 

Proof. Write the complement of the range of G as Ujli, where are disjoint and either of the 
form [ai,bi) or (ai,bi). If x G (0,1) is in the range of G, then G{G~^ (x)) = x, otherwise, if 
x G li, G{G~^ (x)) = bi. By (a), 6j ~ Oj if Oj 0. The last sentence in the statement is then 
proved, and the first follows. □ 

Lemma 3.2. With c\ as in Lemma 3.1, for rj < 1 and j > 2, 

P{G{q,)>r,G{qj))<{l-c,riy-\ 

Proof. By Lemma 3.1 and remarks preceding it, 

P ( > r]G{qj) ) < P (ti > c^rjtj ) = {1 - c^n f'^ . □ 

Lemma 3.3. lim^^oo P {qi < G~^{s/n)) = 1 — e~*. 

Proof. Fix an £ > 0. First, by monotonicity of G~^, ti < s/n implies qi < G~^{s/n). Second, 
by Lemma 3.1 and the monotonicity of G we have that, for large enough n, qi < G~^{s/n) 
implies ti < G{G~^ (ti)) = G{qi) < G{G~^ (s/n) < (1 + e)s/n. These give the inequalities 
P{qi < G-\s/n)) > 1 - (1 - s/n)", and P{qi < G'^s/n)) < 1 - (1 - (1 + e)s/n)". The 
statement of the lemma now follows upon first letting n ^ oo and then £ ^ 0. □ 
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Remark. It follows from Lemma 3.3, and the fact that G{x) = o{x^) near x = 0, that n^^^qi —>■ 
oo as n — oo. 

Lemma 3.4. With high probability qi/q2 is bounded away from 1 as n — ^ oo. More precisely, 
for every rj > there is a S > such that P{qi < {1 — 5) q-z) > 1 — ?7 for large enough n. 

Proof. It follows from Lemma 3.1 that for every 77 > there exists a (5i > so that the following 
implication holds for t2 < Si: if G{qi) > (1 — di)G{q2)) then ti > {1 — r])t2- Furthermore, 
by the assumption (a), there exists & S G (0, (5i) so that, for t2 < (5, gi > (1 — 5)q2 implies 
G{qi) > (1 - 6i)G{q2). Therefore, 

P(gi > (1 - 5)q2) < > (1 - V)t2) + P{t2 >S)=ri + P{t2 > S), 

and the proof is concluded since t2 —>■ a.s. □ 

Lemma 3.5. Yli ~^ ^ as n 00. 

Proof. For any fixed k we have J2'j=i l^/lj — k/nqf 0. Also, J2^=k+i < + l 

a.s. for large n. 

Let (5 > be given. By the above paragraph, it suffices to show that 

lim sup P ( '^^ ^ ^] 

will be arbitrarily small for sufficiently large k. Now, from the assumption (b), it follows that 
for some 77 > wc have > riG{qk+i) whenever qi > 5qk+i and qi < rj. With this rj (which 

we may assume is less than 1) we have, from Lemma 3.2, 

P {-^ X^) < (1 - cirif + P{qi > ri), 

which is clearly enough. □ 

From now on {<^„} will denote a sequence of random variables satisfying = o{qi). Since 
gi S> n~^/^ we shall assume when convenient that also (pn S> n~^/^. In the statement of the 
next lemma, the expression 0{(pn) could have been replaced by the less awkward o{qi). The 
reasons for the present statement arc that the substitute for this lemma (Lemma 6.2) when we 
consider almost sure convergence will have this form, and that the same sequence will 
appear in later lemmas. 



11 



Lemma 3.6. Let {vn} be a sequence of points in a disc with diameter the real interval [— ^ 
0(</^„), C]. Then 

ra— >oo n (1 + r,' 



J=2 



Proof. Write Vn = {bn — ^)/bn- Then if we recall that ^ = (6 — and pj = b — qj we see that 
b — bn lies in a disc with diameter [0, q\ + 0{tpn)\ and that 



n (1 + rjVnY ™ 



hi{b-q,){l-h + q^) 



If we subtract from this the same expression with 6„ replaced by 6, that is, 

n 



(3.1) 

we obtain 
(3.2) 



1 " 

-E 

J=2 



1 " 

-5;(6-g,)(l-6 + g,: 



J=2 



{bn-b + qjY q 



We shall show that this is o(l). Assuming this for the moment, we can finish the proof by first 
noting that we may, with error o(l), start the sum in (3.1) at n = 1 since qi S> n~^/^, and then 
(3.1) has the a.s. limit 

'b\b-q)il-b + qy 

{l + rO' 



It remains to show that (3.2) is o(l). If we replace the numerator 6^ on the right by 6^, the 
error is o(l), since Yl QJ^ bounded. If we make this replacement then what we obtain 

is bounded by a constant times 



n 



E 

J=2 



{bn - bf - 2{bn - b)qj 



qKbn-b + qj) 



A2 



Since \b — 6„| < qi + 0{ipn) = Qi + o{qi) it follows from Lemma 3.4 that |6„ — 6 + qj] is at least 
a constant times qj for large n and so the above is at most a constant times 

1 l&n - b\ ^ 2^ 9l 



J=2 



and by Lemma 3.5 this is o(l). □ 



We denote 
(3.3) 



1 — a/a'^. 



(1 - b) a 

b^e 



1/2 
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Lemma 3.7. We have Un = —1^1^ + Pn~^^^ + o(n~^^^) as n —>■ 00. 



Proof. We show first that u„ > ^ cannot occur for arbitrarily large n. If it did, then we would 
have, using equation (1.1) for n„, 

(C-l)2 - K-l)2 - n (l + nO' ^(l + r.O^' 

It follows from the remark following Lemma 3.3 that the first term on the right is o(l) and from 
Lemma 3.6 that the second term on the right has limit 

since wc are in the composite regime. This contradiction shows that < ^ for sufficiently large 
n, and so Un G [—1^1^, By Lemma 3.6 again, 

n ^ (1 + rjuy " («-l)2 ^"\(l+r^)2 
Therefore the equation (1.1) for becomes 

n(l + rin„)2 U - 1) + / 

Since ri = 6/(1 — 6) + o(l) we find that the solution is as stated. □ 




Next, we see how c„ behaves. 



Lemma 3.8. We have Cn = c(a, F) — 6 qi + o{qi) as oo, where 9 is given in (3.3). 



Proof. Write 
(3.4) 



1 a X ^ fjUn oc 

l — Un n 1 + rjUn n 1 



a riUj, 



By Lemma 3.7, the last term above is 0{n ^/^). Equation (1.1) tells us that 



+ riUr, 



du 



I g rju \ 

\l — u n 1 + rju J 



0, 
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and so 



d I I a TiU 



du \ 1 — u n 1 + rju j u=u„ n (1 + riu„)2 ri/?^ 6/3^ 



By Lemma 3.6 and its proof, with an error o(l) the derivative of the expression in the 
parentheses above equals in ["Un,^] what it equals at u = ^, so the above holds with u„ replaced 
by any point in this interval. Prom this and (3.4) we get 



We have 



where we have used the fact that qi 2> ri~^/^. Thus 

a(l -h) , , 

Cn = C(0 53^51 + 0(gi). 

Finally, as ((6 — p)^) < oo, we can use the central limit theorem to conclude that c(^) 
c{a,F) + 0(n~-^/^), which completes the proof. □ 



Remark. Lemmas 3.3 and 3.8 show that Theorem 2 follows from the part of Theorem 1 on 
convergence in probability. 

4. Preliminary lemmas II: steepest descent curves. 

Now we go to our integrals (2.1) and (2.2). We are not going to apply steepest descent with 
as the main integrand, but rather with the function ■i/^i which is with the factor 1 + r-\_z 
removed. It is convenient to introduce the notation 

n 

M^,c) = n(l + rjz) {z - ir z-^'-^^^, 

where c > 0. (This parameter is not to be confused with the time constant c = c(a, F) defined 
earlier.) Thus ipi{z) = ipi{z,Cn) in this notation. We also define the integrals 

/+(c) = 2^/ {l + nz)Mz,c)dz, r{c) = ^-j (l + riz)-Vi(2;,c)-^z-'dz. 
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(Since I~^{c) = when c > 1 we always assume that c < 1.) Notice that these are exactly the 
integrals (2.1) and (2.2) when we set 

c = Cn + {h' + j + k)/m. 

Since j,k >0 and we will eventually set h' = sn^^^, we may also assume that 

(4.1) c>cn-0{n-'/^). 



To apply steepest descent to I^{c) we must locate the critical points and determine the 
critical values of ijji{z,c). Thus we define 

ai{z,c) = — log tlJi{z,c), 
m 

so that 

c- 1 



a[{z,c) = -y — + - + 

^ n ^ 1 + r,-2; z-1 



As before, if the parameter c does not appear we take it to be c„, e.g., C7i{z) = ai{z,Cn)- So 

a[{z) = — log t/ji{z) = a'{z) - -7-^ — • 
m n 1 + riz 

Using a'{un) = a"{un) = we get from the above and Lemma 3.7 that 

(4.2) ^;(^„) = __^(l + o(l)), ^;'(^„) = ^(i + o(i)). 



To determine the critical values of ai {z, c) let us first find the value of c for which its derivative 
has a double zero. (This is the analogue of the quantity for a{z).) For this we use the analogue 
of (1.1) and (1.2) but where the terms corresponding to j = 1 are dropped from the sums. If we 
call the solution of (1.1) u and set c = c(u) then a[{z, c) has a double zero at u. In analogy with 
Un, we know that u is to the right of and within 0(n~^/^) of — r^^. As for c, we use Lemma 3.8, 
its analogue where the sums in (1.1) and (1.2) start with j = 2, as well as Lemma 3.4, to see 
that to a first approximation 

c = c„ - e{q2 - qi) 
and that g2 — 9i S> n~^/^. Prom this and (4.1) we see that c > c. 



Using subscripts for derivatives now, we have 
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and we want to see how the critical points of (Ti{z, c) move away from u as c increases from 
c. (Here we take u~ < The function aiz{z,c) vanishes at u and is otherwise positive in 

(— r^^, 0). It follows that for c close to but larger than c we have u~ < u < u^. Differentiating 
o'lz {uf ) c) = with respect to c gives 

(4.3) = (7u^ {uf , c) ^ + {uf , c) = {uf,c)^ + ^. 

ac ac Uc 

Since < it follows that du^ /dc ^ 0, and so each of is either a decreasing or increasing 
function of c for c > c. Prom their behavior that we already know for c close to c we deduce 
that u+ increases and u~ decreases as c increases. In particular, u~ is even closer to — than 
u. 

We remark that from (4.3) and the signs of du^ /dc we deduce 

(4.4) (Tizz{ut^c)>{), ai^^ {u- , c) < 0. 

Next we shall determine the asymptotics of the critical values a{u^,c). The sequence {<^n} 
is as described before Lemma 3.6. 

Lemma 4.1. For c — Cn = 0{ipn) 

(4.5) a,{ut,c) = a,i-r^\ c) - (^c - c,, + ^{1 + o{l))n-'/'^ ' 
and for all c> Cn 

(4.6) ai{u+, c) < ai{-r^\ c) - jyn-^/^ _ + 0(„-i). 
for some rj > 0. Moreover for all c 

ai{u-,c) > ai{-r^^,c) + (pl 

when n is sufficiently large. 

Remark. In these and analogous inequalities below we think of ai as actually meaning 5R(7i. 
Proof. Consider first the case c = c„. We have 

<Tl{Un + C) = <^l{Un) + <T'i{Un) C + C f (1 - t) (Ti (t^n + ^C) dt. 

Jo 
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If C = 0{ipn) then it follows from Lemma 3.6 that a1{un + tQ = (t" + o(l). Hence, by (4.2), 
we have for such ( 



(4.7) a,iu^ + = ^i(^n) - +[2^2+ ) C • 
This has zero derivative for 

C=^(l+o(l)) 

and it follows that 

(4.8) ul =un + ^il + o(l)) = -rr^ + + 

(This critical value must be 'u+ rather than since the latter is within 0(ri~^/^) of — r^^.) 
Prom this and (4.7), taking C = -r^^-u^ = -(/3+o(l))n"^/2 andC = u+^-Un = (/3+o(l))n"^/2 
and subtracting, it follows that 

(4.9) ai«) = ai(-rr^) - 2(a + o(l))n-^ 



To determine the behavior of u'^ and ai (u J , c) for more general c we assume first that 
c = c„ + o(l), ii+ = u„ + 0(¥?„) = -rf ^ + 0{(pn)- 

Then 

fTlzz(u+, c) = a"{un) ^ " ^ + ^^^^ 

by (4.2). Therefore (4.3) gives 



= -(/3V« + o(l))K = r-A(l + o(l) 
etc q: 

whence 



/32 

«c =<+'^17-(c-Cn)(l + o(l)) 



+ ^(1 + o(l)) + ri^(c - c„)(l + o(l)), 
n a 



by (4.8). This holds if c— c„ = 0{ipn) since this assures that n+ = n„ + 0((^„). The above gives 
(4.11) log(-tx+) = log(-rr^) - 2r^P{l + o{l))n-'/^ - r?^(c - c„)(l + o(l)). 



a 

(Again, real parts are tacitly meant.) 
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To determine, ai (n^ , c) we use aiz {u'^ , c) = to deduce 
(4.12) ^(7i(u+,c) =logn+. 

Wc continue to assume that c — c„ = 0{fn) so our estimates hold. Integrating (4.12) using the 
first part of (4.10) gives (since u'^^ — > — r^^) 

1 /S^ 

C7i(m+,c) = 0-1 (U+) + (C- C„) logU+ - 2''l— (C- Cn)^ (1 + o(l)) 

= f^i(-r-r') - 2(a + o(l)) + log(-rr')(c - c„) 

- 2ri^(c - c) n-i/2(l + o(l)) - ^r?^ (c - c^)^ (1 + o(l)), 

2 a 

by (4.9) and (4.11). This gives (4.5). 

For all c > c„ we use the fact that log(— n^) is a decreasing function of c, since u'^ increases, 
and integrate (4.12) with respect to c from c„ to c, which gives 

(^i{u^,c) < cri(n+ ) +log(-n+ )(c- c„). 

Using (4.9) and (4.8) give (4.6). 

For the lower bound for ai{u~,c), we assume first that c < c„. By (4.1) this implies in 
particular that c — c„ = 0(n~^/^). Now cri(z) is decreasing on the interval {u~ , «+) and 
n+ — -u" S> tpn- To see the last inequality, note that, from Lemma 3.6, aizz{un + CjC) 7^ for 
C = 0{ipn) and c — c„ = o(l). Therefore aiz{un + (,c) can vanish for at most one such ( and, 
since — Un = 0{ipn), we must have Un — u~ ^ (pn- 

Take any sequence ipn = o{qi) and write 

ai{u~,c) > ai{u+ - (pn,C) = ai{u+ - (pn) + {c - Cn) log{ipn - li+). 

(As usual, we imagine real parts having been taken.) If we apply (4.7) with ^ = n+ — n„ and 
with C, = — — Un and subtract, we obtain 

U{ut - ^n) - 0{ut) = ^n-l/Vn(l + o{l)) + ^ (-2^„(«+ - «„) + )) (1 + ^(1)). 

By subtracting the first parts of (4.10) and (4.8) we see that this equals 
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Since S> n 



-1/2 



, as we may assume, we obtain 



for some 77 > 0. Also, since c — c„ > —rin for some 77 and log(l — ipn/u^) is positive and 
0{fn) we have 

(c - c„) log((/?„ - n+) > (c - c„) log(-u+) - rin~'^/'^ipn- 



Putting these together gives 



for some 77 > 0. 

This was for c < c„. For c > c„ we use what we get from (4.12) by replacing + with ~, 
subtracting the two, and integrating. Together with using the already proved inequality for 
c = c„ this gives 



CTi (n^ ,c)-ai (n+ ,c)>ri(pl+ / \og{u^ / u+ ) dc. 
The logarithm is nonnegative. Hence ai{u~ ,c) — a"i(?x+,c) > rj^p^ for all c. 

If c — c„ = 0{ipn) then using this and (4.5) give 

ai{u~ ,c) > ai{-r^^) + \og{r^^){c - Cn) + wl- 
with a different rj. If c > c„ we use 

C7i {u- ,c)-ai («- ) = / log(-u^ ) dc. 

Since ?x~ is decreasing and is less than —r^^ when c = c„ this gives 

cri(n-,c) > cri(u- ) +log(rj7^)(c-c„) 

> (71 «) + log(ri-^)(c - c„) + ifl. 

Combining this with (4.5) for c = c„ shows that 

holds for these c as well. Since {(pn} was an arbitrary sequence satisfying y?„ = o(gi) the last 
statement of the lemma follows. □ 

Next we consider the steepest descent curves, which wc denote by C^(c) corresponding to 
the integrals /^(c). It follows from (4.4) that C~^{c) passes through tt+ because on the curve 




19 



|'0i(2;, c)| has a maximum at that point; similarly, C~(c) passes through u~ . We have enough 
information to evaluate the portions of these integrals taken over the immediate neighborhoods 
of these points, but we also have to show that the integrals over the rest of the curves are 
negligible. This requires not only that the integrands are much smaller there, which they are, 
but also that the curves themselves are not too badly behaved. 

To see what is needed, let r='= be arcs of steepest descent curves for a function p, curves on 
which SJp is constant. In analogy with our C^(c) we assume Sip is increasing on T~ as we move 
away from the critical point and decreasing on r+. If s measures arc length on r='= we have for 

(4.13) ^ = =fJ4^. 

as p'{z) 

If the arc goes from a to 6 then 

\p'{z)\ ds = Tj^ P'{z) dz = T{p{b) - p{a)). 
Hence the length of r='= is at most 

mm^eri \p'[z)\ 

This is to be modified if p' has a simple zero at z = a, for example. In this case we replace p'{z) 
by p'{z)/{z — a). (This is seen by making the variable change z = a + ^/^.) 

Our goal is Lemma 4.5 below. In order to use the length estimate (4.14) to deduce the bounds 
of the lemma, we must first locate regions in which our curves are located, and then find lower 
bounds for a[{z,c) in these regions. (Upper bounds for |(Ti(z,c)| will be easy.) These will be 
established in the next lemmas. 

For r > define n(r) = #{j : rj > r}. 
Lemma 4.2. The curves C^(c) lie in the regions 

f CTl I 

<z:\ arg(r-^ + z)\ < tt — -— \ 

I an{r) + cn ) 

for all r and in \z + r2^\ > 5n~^ if S is small enough. 

Proof. For a point z on either of the curves, say in the upper half-plane, we have 

n 

CTT = — arg(r~^ + z) + aTg(z — 1) + (c — 1) arg z 
" i=2 

^ an(r) , , , _i , 

> ^ arg(r + z) + c arg(r + z), 
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which gives the first statement of the lemma. For the second, observe that if C = 0{(pn) then 
(Ti(r^^ + C, c) = a/n( + 0{l). This shows, first, that hes to the right of the circle \(\ =dn~^ 
if S is small enough and, second, that l/a[{z,c), thought of a vector, points outward from this 
circle if 6 is small enough. Since a point of C~(c) moves in the direction of l/a[{z. c) as it moves 
away from u~ (see (3.7) of [GTW2]), the curve can never pass inside the circle. Therefore the 
entire disc |C| < 6n~^ lies to the left of C~(c). This gives the second statement for C~(c) and 
it follows also for C+(c) since this is to the right of C~(c). □ 

The next lemma, together with (4.13) and the length estimate (4.14), will imply that for z 
large the curves will move in the direction of z and are well-behaved. If we take any f < 5/(1 — 6) 
then a positive proportion of the rj are greater than f and so by Lemma 4.2 the curves lie in a 
region 



for some d > 0. 

Lemma 4.3. We have za'i{z,c) ^ c + a as n —> oo and z — > oo through region (4-15). 
Proof. We have 



and it suffices to show that the last term tends to as n — *■ oo and z ^ oo through region (4.15). 
If z is in this region and r < f/2 then |1 + rz| > (5(1 + r\z\) for another 5. The same bound will 
hold for all r < 6/(1 — 6) if z is large enough. Choose M large and break the sum on the right, 
with its factor n~^, into two parts, the terms where rj\z\ < M and the terms where rj\z\ > M. 
We find that its absolute value is at most 



The first term tends to as 2; — >^ oo while the second could have been arbitrarily small to begin 
with. □ 

Remark. If P{p = 0) is positive then the above has to be modified. We replace c + a by 
c + aP{p > 0). 

Because of the above lemma we need only consider 2; in a bounded set. We use the fact that 
by Lemma 4.2 with r = r2 our curves lie a region 



(4.15) 



{z : I arg(f"^ + z)\ < 7r(l - 5)} 



za[{z,c) =c + a + 0{n~'^) + 0{z~^) + - 




n ^ 1 + ' 

i=2 3 



n-i(n-n(M/|z|)) + -i-. 



(4.16) 



{z : |arg(r2 ^ + z)\ < 7r(l - Sn~^), ^ + z\ > Sn~^} . 
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Lemma 4.4. For all z in any bounded subset of the region (4-16) we have 

{z - U-) {z - u+) 



\a[{z,c)\ > 6n-^ 
for some S > independent of c. 



z{z - 1) 



Proof. To obtain the lower bound we write 

-1 

a sr^ 1 

!, S3,-- - ,Sn; Z) = - > 



(s:z 



(S2, 



n '■ — ' So + z z—\ 

j=2 ^ 



1 c- 1 

+ 7 + ■ 



Of course a[{z,c) = <?i'(?'2 ^ > ''a ^ > - - ' , ^ ) • Think of ^2 = ^2^ ^ ^ fixed, and consider the 
problem of finding inf |(/>(s; z)\ where S3, • • • , s„ are subject to the conditions 



Sj>S2, 4>{s;u^ 



0. 



If we take sequences so that the inf is approached in the limit, then some Sj may tend to infinity, 
others may tend to S2, and the rest, if any, tend to values strictly greater than S2. Thus our inf 
is equal to the minimum of \(j){s; z)\, where (f) now has the form 



n ^ 

, I .a Ua 1 

9\S2-, S3, - - - , Sn'; = - V 7 \ T + 

n ^ s^-V z z — 1 
j=2 ■> 



c- 1 



with n' < n, ^rij = n — 1, and the Sj with j > 2 satisfying Sj > S2 and the constraints 
0(s;u±) = O. 

Notice that the minimum cannot be zero since (/)(s; z), thought of for the moment as a function 
of z, has n' finite zeros. It has zeros at and one between each pair of consecutive —Sj since 
all the coefficients of l/{sj + z) are positive. This accounts for all n' zeros, so our z cannot be 
one of them. 

We apply Lagrange multipliers to find the minimum oi \(f){s; z)\'^ over S3, • • • , s„', achieved at 
interior points. There are two constraints, hence two multipliers A and /x. If p + iq is the value 
(f){s; z) where its absolute value achieves its minimum, then the equations we get are 



^{p — iq) 



A 



+ 



{sj + zy (sj+uc)^ {sj+uty 



where we have divided by the factor nj appearing in all terms. This is the same sixth degree 
polynomial equation for all the sj. It follows that there are at most six different Sj. Assuming 
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there are exactly six (if there are fewer the argument is the same and the final estimate is better) 
we change notation again and write these as S3 , • • • , sg so that the minimum is achieved for 



8 1 1 

a \ rii 1 c — 1 



(P{s2, S3,-- - ,ss; z) = - y — f \ + 



n Sj + Z 

j=2 ■> 

with other rtj . 

This has eight zeros. Two of them are and the other six, lying between consecutive — Sj, 
we denote by lii, • • • , tig. We have the factorization 

Sis- z\ = 1 ~ g - ^j) - ^t) nLi(i - ^M) 
u-,ut z{z-\) n-=2(i-^/^i)' 

and it remains to find a lower bound for this. Near z = we have o'^iz, c) = (1 — c)z~^ — 1 + 
a (r) + 0(1), so if c is close to 1 then (1 — c)/ti+ = 1 — q (r) + o(l). In particular this is bounded 
away from zero. Thus the first factor above is bounded away from zero. As for the factors in the 
products, observe first that each factor \ — zjsj is bounded since z and all factors 1/sj are. For 
the others, we use again the fact that the curves lie in a region (4.16). In any bounded subset 
of this region each |1 — zlui\ > r)n~^ for some t] > 0. (If z is in a neighborhood of this is clear 
since each Ui < 0. Otherwise write 1 — z/ui = z{z~^ — u~^).) Therefore the product of these is 
bounded below by a constant times n~^. This completes the proof. □ 

Now we can show that the curves C^(c) are not too badly behaved. 

Lemma 4.5. For some constant A> the length of C~^{c) is 0{n^) and 

[ \z\-^ \dz\ = O(n^). 

Jc-{c) 



Proof. It follows from Lemma 4.3 that C~^{c) lies in a bounded set. For, this lemma implies 
that the vectors l/a[{z,c) point outward from a large circle \z\ = R, and since by (4.13) C~^{c) 
goes in the direction opposite to l/a[{z,c), a point of the curve starting at can never pass 
outside the circle. Also, some disc \z\ < S{1 — c) is disjoint from C~^{c) because l/a[{z, c) points 
outward from a small enough circle \z\ = 6{1 — c) and so C~^{c) cannot cross into it. It follows 
that cr[{z, c), and so also (Ji{z, c), is bounded on any portion of C^{c) close to z = 0. A similar 
argument shows that some disc |;z — 1| < S lies entirely inside C^{c). Finally, we know that u~ 
is within 0{n~^^'^) of — and if C = o((/i) then o"i(r^^ + (,c) = a/nQ + 0(1)- In particular 
u~ lies in a region \C,\ > 5n~^ for some 5 > 0. Since also a'l = —a/nC,"^ + 0(1), by Lemma 3.6, 
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we deduce that ai{z^c) = 0{n) when \z — u~\ < dn~^/2, thus for such z we have ai{z,c) = 
(7i(n~,c) + 0{n\z — But it follows from Lemma 4.1 that ai{u~,c) — ai{u^,c) > cp^, 

and then, since = o((^^), ai{u^,c) < ai{z,c) for \z — u~\ < Sn~^/2. As the maximum of 
ai{z, c) on C~^{c) occurs at n+, this shows that the distance from C~^{c) to u~is at least Sn~^/2. 
With these facts established we use the lower bound of Lemma 4.4, the length estimate (4.14) 
(extended as in the remark following it), and the obvious upper bound for |c7i(2;, c)| in the region 
(4.16) to deduce that the length of C+(c) is 0(n"^) for some constant A. 

As for the integral over C~(c), we observe that, since c < 1 and cm is an integer, 1 — c 
is at least a constant times n~^. Since C~(c) lies outside a disc \z\ < S{1 — c), we have 
z~^ = 0{n) on C~(c). A lower bound for the distance from C~(c) to is obtained using the 
fact that ai{u~,c) — ai{u^,c) > (p"^. Since a[ is bounded in a neighborhood of u'^, we have 
ai{u~,c) > ai{z,c) for Iz — u+l less than ip'^ times a sufficiently small constant. This shows 
that C~(c) is at least this far from u^. We apply the other bounds as before; we think of the 
integral over the portion of C~(c) outside a large circle as the sum of integrals over the arcs 
from Ofc to Ofe+i where is the point of C~(c) where \z\ = k. Lemma 4.3 and (4.14) are used 
again here. □ 

5. Asymptotic evaluation of the integrals. 

We evaluate I~^(c) first when c — c„ = 0((/?„). Then aizziu^jc) = a//3^ + o(l) and so if we 
set z = + ( we have 



a 



2 



ai(z,c) =(7i(u+,c) + 2^(l + o(l))C 

as long as C = 0(<^n)- If ICI = <Pn then the real part of the second term above is less than a 
negative constant times ip'^ and, since this real part decreases as we go out C"'"(c), it is at least 
this negative whenever \(\ > ip^. If we recall that this gets multiplied by m in the exponent 
and the fact that C^{c) has the length at most a power of n (by Lemma 4.5), wc sec that the 
contribution of this part of the integral is O ^e"*'^^"^ . It follows from Lemma 3.3 
and assumption (c) that with high probability qi S> logn/n^^'^ , and we could have chosen ip^ 
to satisfy this also. Thus, with error o(e™''^("<^''^)) the integral I~^{c) is equal to 



- / (l + ri(t^+ + C))e("/2^')(i+°«)^'dze--^(< 



■ c) 



(since am = n). Since ipn J?> n"-"^/^, in the limit after making the variable change C, n~^/'^C, 
the integration can be taken over {—ico, ioo) (downward really, but we can reverse the directions 
of integrations), the linear factor contributes zero, and by (4.10) 

l + riu+ = ri (2(3n-^/^ + ^ic-Cn)+o{n-^/^ + \c-Cn\)^ . 
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Thus the integral is asymptotically equal to /3\/2Trin times the above and, by (4.5), 

7+(c) = ^n-i ( 2 + ^ni/2(c - c„) + o(l + n'/'\c - c^A 
V27r V " / 



X 



This assumed that c — c„ = 0{ipn)- For all c > c„ we use the second part of Lemma 4.1 and 
again the fact that C+(c) has the length at most a power of n. We deduce 

/+(c) = O (^V'l(-rr\c) e"''"'^' (c-c„)+0(logn)^ 

for c > c„. 

For the integral over C~(c) we use the last part of Lemma 4.1 and the second part of 
Lemma 4.5. These imply that the integral over C~ is 



O 



(V'i(-rr\c)-ie-"'^"+o('°s")) =o(V^i(-rr\c)). 



But our integral for I~ (c) is not taken over C~ (c) . Recall that the original contour must 
have all the —rj^ on the outside whereas —r^'^ is inside (more precisely, on the other side of) 
C~{c). Therefore if we deform the contour to C~(c) we pass through the pole at —r^^- Thus 

/-(c) = n V'i(-rj"\c)"^ + o{ipi{-r^^,c)). 



Now recall that in I~^{c) we set c— c„ = h'+j+£, in I~{c) wc set c—Cn = h'+l+k and then we 
sum over f, to get the matrix product. Recall also that 'il}i{—r^^,c) = ipi{—r^^) (— ri)~"*^°~'^"\ 
The factors (— n)""**^^"^") in /"'"(c) and (— ri)"*(^~°") in /~(c) will combine to give (— n)™-*^^"-^) 
which can be eliminated without affecting the determinant. It follows that we can modify the 
expressions for /^(c) by removing these factors. We can also remove the factors ipi{—ri^)^^ 
since they cancel upon multiplying. Thus our replacements are 

V27r V " / 
if c- c„ = 0{(pn), and 

if c > c„. Furthermore, I~{c) ri + o(l). 
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Recall next that we set h' = sn}/'^ and in /"""(c), c = c„ + sv}^"^ + [xn-*^/^] + [zn^/^J, so that 

c — Cn = {s + x + z + o(l))n^/^/m = Q!(s + X + z + o(l))n~^/^, 

and eventually we multiply by n because of the scaling. Take first the case c — c„ = O ((/?„), that 
is, a; + z = 0(n^/^(^„). Since m = n/a and rif3 = (1 + o(l)) the modified I^{c) equals 

4Sn-l(2T + S + X + Z + o{l+X + y))e-^^r+s+x+z+o{l)f/2r\ 



On the other hand, / (c) is equal to ri with error o(l). The result of multiplying these together, 
multiplying by n, and integrating with respect to z over (0, oo), is asymptotically equal to 

(5.1) ^_g-(2r+.+x)V2r^_ 



This holds for c — c„ = O (</?„). If c — c„ > we have, for our modified /"'"(c), the estimate 

Integrating the square of this over a region x + z = 0(n^/^) will give o(l). 

It follows that the matrix product scales to the operator on (0, cxd) with kernel (5.1). This is 
a rank one kernel so its Fredholm determinant equals one minus its trace, which equals 

1 f ^-xV2r\ 



'2irT J -oo 

This establishes the convergence in probability statement of Theorem 1. 

Remark. One could rightly object that to scale a product to a trace class operator we should 
know that each factor scales in Hilbert-Schmidt norm. In our case the second limiting kernel is a 
constant and the product is not even Hilbert-Schmidt. But we could have multiplied the kernel 
of the first operator by (1+,t) (1 + -z) and the kernel of the second operator by {\ + z)~^ {l+y)~^ . 
This would not have affected the determinant of the product, both operators would have scaled 
in Hilbert-Schmidt norm and the product would have scaled in trace norm to the rank one kernel 

1 g-(2r+s+x)V2r2 1 + 3^ 



which has the same Fredholm determinant. 
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6. Almost sure convergence. 

What is needed, and all that is needed, is an "almost sure" substitute for Lemma 3.6 under 
assumptions (a') and (b'). We begin with a lemma on extreme order statistics of uniform random 
variables, part or all of which may well be in the literature. 

Lemma 6.1 Let a> 1 be arbitrary. Then, almost surely, 

t.>-^, ^<1 ' 



n log n t2 log n 

for sufficiently large n. Here, r) is a positive constant depending on a. 

Proof. We use the notation tnj for our tj to display their dependence on n. We have 

P{tn,i < 5) = 1 - (1 - 5)" ~ n(5 if nS = o(l). 

In particular 

2-'=\ 1 

P U2M < 



It follows that, a.s. for sufficiently large k we have 



t2k I > 



Take any n and let k be such that 2^^ ^ < n < 2'^. Prom the above we have, a.s. for sufficiently 
large n 

tn,l > ^2*^,1 > -J-T- > — i 

fe" nlog n 

for some 77. 

For the ratio we use the fact that 

(6.1) P I - s] =1-{1- 6y ^ j6 if jS = o(l). 

Now suppose that 

(6.2) ^>1--^ 

tn.2 log"n 

and let k be such that 2*^"^ < n < 2*^. Take any J (which will eventually be of order log A;). 
Then there are two possibilities: 
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(1) t2kj < tn,i for all j < J; 

(2) t2kj > tn,i for some j < J. 

Consider possibility (1) first. Let G„ be the event that tn,i < aloglogn/n. By Ex. 4.3.2 of 
[Gal], P{Gn eventually) = 1. Moreover, 

P{{t2kj < tn,i for all j <J}n Gn) < -P(^2^J• < 21oglogn/n for all j < J) 



for some constant A. If J = Blogk then the bound above equals e~^^^°^^~^^^°^'' , so if we 
choose B large enough the sum over k of these probabilities will be finite. With this J, (1) can 
therefore a.s. occur for only finitely many k. 

Next consider possibility (2) and let j be the smallest integer < J such that t2kj > tn,i- 
Then t2kj < t„,2 and = t2k^£ for some £ < j. It follows that t2fc,j-i/i2'=,j > in,i/tn,2 and 
by (6.2) this is at least 1 — C/Zc", for some constant C (which will change from appearance to 
appearance). Therefore, by (6.1), 



It follows that (2) and (6.2) can happen together only for finitely many n. The upshot is that 
a.s. the inequality (6.2) can occur for only finitely many n, which completes the proof. □ 

We are now ready to prove our substitute for Lemma 3.6. Recall that we can set Qj = G-^{tj). 
The assumption (a') implies that G is continuous near 0, so that G{G~^ (x)) = x for small x. 



Proof. Prom the proof of Lemma 3.6 we see that we want to show that, for some sequence (/p. 
as described, we have a.s. 




P((6.2) and (2) both happen) 

< ^(^2^i-lA2^J• > l - C/k'' for some j < J) < CJ^/k"" < Clog^ k/k". 



Lemma 6.2. Suppose (a') and (h' ) are satisfied. Then there exists a sequence ipn 3> log n/n-^/^ 
such that a.s. for any sequence {vn} lying in the disc with diameter the real interval [—r^^ — 
0{ipn)^ C] we have 





n 
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Assumption (a') implies that 



y - \G-^y) 

when X < y are small enough. Therefore, it follows from the second part of Lemma 6.1, that 
a.s. for large n, 

(6.3) ^ < 1 ^ 



for another constant rj > 0. Set 



q2 log n 



I 1 r] 
Vn = t;-, — —Q2- 
2 log n 



Let us show that ^„ ^ logn/n^^'^ . Assumption (a') implies that G~^{x) is at most a 
constant times x^^"', thus the fact that ti = 0(log log n/n) shows that qi is at most a constant 
times (log log n/n) ^Z'''. Furthermore, assumption (b') gives, with a slightly smaller v, .x^ S> 
G{x) log'^ x~^. Applying this with x = qi = G~^{ti) and using the first part of Lemma 6.1 gives 

nlog n 

We therefore deduce that 

(6.4) -log^'-'^n 

n 

for a slightly smaller v than in (b'). By (6.3), the same holds for q2 and so 

V'^»-log''-=^»n 
n 

and tpn ^ log n/n^/^ as long as — 3a > 2. Since a > 1 is arbitrary the requirement becomes 
v > 5. But from (a') and (b') we see that necessarily 7 > 2, so that v > 8. 

If j > 2, then (6.3) and the inequality q2 < Qj imply that 

g,-(5i + V.)>^l^?,. 



We take for any sequence satisfying 

logn 

—r-pr <^(Pn^ Wn- 

At this point we follow the proof of Lemma 3.6 to see that the expression 

1 2a n 

(6.5) l^y-^ 
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needs to go to a. s. to conclude the proof of this lemma. This is what we will demonstrate. 

For any kn, if wc separate the sum in (6.5) over j < kn from the sum over j > kn, we see 
that (6.5) is at most 

(6.6) 2-fc„ + log n 2^^- 



Wc first determine kn so the second term in (6.6) goes a.s. to 0. By strong law, n ^^Qj ^ 
{q~'^) a.s., so log^" ^Qi/lkn^-i needs to go to 0. We have, for each 6 > 

P f log- > P f > ^) 

- - Vlog-ny J V Vlog-ny J 

This is summable over n if we choose 

A;„= Llog«n(log2«n)^J + l. 

With this choice, the second summand in (6.6) therefore goes to a.s. 

On the other hand, the first term in (6.6) is with the same choice of fe„ at most a constant 
times 

log(27+3)a^ 

and from (6.4) this is o(l) times log'-^'*'"'"^^""'' n. Since a > 1 was arbitrary and u > 2'y + 4, we 
can make (27 + 4)a — u < and then the first summand in (6.6) goes to a.s. This completes 
the proof. □ 



With this lemma in place of Lemma 3.6 the reader will find that all subsequent limits and 
estimates in Sections 4 and 5 will hold almost surely, thus giving the second statement of the 
theorem. The reason our sequence had to satisfy (pn ^ logn/n^^'^ is that errors of the form 
O (e~"''^""'"'^('°S"') J appeared in the evaluation of /^(c) and these had to be o(l). 
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